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Abstract —  In  this  paper  a  wavelet  based  method  for  SAR  data 
denoising  and  compression  is  presented.  An  unsupervised 
stochastic  model  based  approach  to  image  denoising  is  presented. 
SAR  image  is  modeled  in  wavelet  domain  Gauss  Markov 
Random  field  and  noise  is  considered  as  Gaussian  with  unknown 
variance.  The  parameters  are  estimated  from  incomplete  data 
using  mixtures  of  wavelet  coefficients,  and  expectation 
maximization  algorithm.  The  expectation  maximization 
algorithm  is  used  to  efficiently  compute  a  maximum  a  posteriori 
estimate.  Observed  wavelet  coefficient  is  estimated  using  inter 
and  intra  scale  of  wavelet  coefficients  to  estimate  image  and  noise 
model  parameters.  Presented  wavelet  based  method  efficiently 
removes  noise  from  SAR  images.  The  second  step  is  to  design  an 
entropy  coder  that  efficiently  codes  despeckled  image.  The 
texture  parameters  obtained  at  the  despeckling  stage  are  used  in 
the  compression  process.  The  image  coder  is  tested  on  X-SAR 
data  with  and  achieves  comparable  compression  results  with  the 
wavelet  based  state-of-the  art  coders  for  SAR  data  compression. 
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I.  Introduction 

Over  the  past  few  years  all  kinds  of  data  presentation  has 
been  intensively  investigated.  The  SAR  (synthetic  aperture 
radar)  data  provided  by  X-SAR  and  RADAR-SAT  missions 
have  10-20  meter  resolution.  SAR  radars  with  1  meter 
resolution  have  been  intensively  investigated.  The  Terra-SAR 
radar,  will  provide  images  with  0.7  meter  resolution. 

In  this  paper  we  will  introduce  a  wavelet  based  denoising 
for  SAR  data.  Many  different  wavelet  based  techniques  have 
appeared  in  recent  times.  It  is  well  known,  that  SAR  data  are 
corrupted  by  multiplicative  noise  called  speckles.  If  the 
speckles  are  removed  from  the  SAR  image,  higher 
compression  ratios  for  data  compression  can  be  achieved  [1]. 
Many  authors  have  shown  that  the  speckles  can  be  efficiently 
removed  in  the  wavelet  domain.  Two  different  approaches  are 
usually  applied  for  speckle  removal;  the  first  one  uses 
maximum  a  posteriori  (MAP)  filtering  of  the  image,  where  the 
image  and  speckles  are  separately  modeled  [2]-[3],  and  the 
second  approach  is  to  use  MMSE  estimator  as  shown  in  [4]- 
[5].  With  other  words,  all  authors  have  introduced  a  shrinkage 
factor  that  predicts  the  coefficient  without  noise.  If  the  image 
model  is  well  chosen,  the  despeckling  algorithm  can  provide 
side  information  about  the  texture  or  even  the  image  can  be 
segmented  into  different  areas.  When  the  image  is  modeled  as 
a  Gauss  Markov  Random  field,  texture  parameters  can  be 
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derived.  Many  algorithms  for  segmentation  have  been 
proposed.  In  [6]  a  Markov  random  field  is  used  for  image 
model  and  expectation  maximization  (EM)  [9]  is  used  to 
estimate  parameters.  Bayesian  tree  structured  image  modeling 
using  wavelet  domain  hidden  Markov  models  is  proposed  in 
[7].  Multiscale  Bayesian  segmentation  using  a  trainable  context 
model  is  proposed  in  [8].  In  [7]  and  [8]  a  multi-scale  wavelet 
trees  are  used  at  modeling  stage.  In  this  paper  we  introduce 
wavelet  based  algorithm  for  SAR  data  despeckling.  The  image 
is  modeled  as  a  Gauss-Markov  random  field  and  a  noise  is 
considered  as  a  Gaussian  noise  with  unknown  variance  and 
mean.  A  Bayesian  approach  is  used  to  determine  a  posteriori 
and  with  expectation-maximization  parameters  are  estimated. 

II.  MODEL  BASED  DENOISING 

Let  y  =  (x,  z)  be  a  set  of  random  data  of  interest,  where  v  is 
observed  data  and  z  is  a  part  of  y  that  is  not  observable. 
Observable  data  x  is  usually  called  incomplete  data  and  y  is 
called  complete  data.  Let  p[x  \  O)  be  the  probability  density 
function  (pdf)  of  incomplete  data  x,  where  O  is  a  set  of 
parameters  that  characterize  pdf  and  p(^y  \  O)  pdf  of  complete 
data.  The  incomplete  data  problem  is  to  estimate  O  based  only 
on  the  observations  presented  by  incomplete  data  x.  The  pdf  of 
incomplete  data  is  then 

piA^)=L,piym  (1) 

An  iterative  maximum  likelihood  (ML)  approach  is  to  find 
estimate  : 

=  arg  max  /7(x  |  O)  (2) 

o 

The  most  general  and  effective  algorithm  for  solving  this 
problem  appears  to  be  EM  algorithm  [9].  The  EM  is  iterative 
algorithm  and  requires  expectation  and  maximization  steps.  In 
the  E-step  of  the  EM  algorithm  we  have  to  compute  the 
functional  for  the  subsequent  maximization  in  the 

M-step. 

Q(<I)  1 )=E  [log  p{y  I  O)  I  X,  O'") ]  (3) 

In  the  M-step  a  maximization  of  functional  is  computed: 
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Q(<i)('>+'))=argmax 0(0 1  O*'’*)  (4) 

<5 

To  apply  EM  algorithm,  we  start  at  initial  <^0)  and  iterate 
algorithm  until  algorithm  converges:  -O'’!  <  £■,  for 

appropriately  chosen  vector  norm  and  value  £  . 

A.  Model  selection 

In  this  paper  the  image  is  modeled  using  Gauss-Markov 
random  field: 


p(w^\Xi,9)  = 


rCXp 


^ 

2 


Ha: 


V  J  J 


(5) 


Speckle  model  can  be  chosen  as  a  Gaussian  process: 


o(h’,  I  X,)  =  ,  ^  exp(-(H’, -x,Y  / 2cr/') 

72^  V  ,  .  I 


(6) 


In  equations  (5)  and  (6)  variances  cr^  and  are  unknown. 
The  parameter  0  that  provides  texture  information  can  be 
estimated  on  the  various  ways.  In  [2]  maximum  a  posteriori  is 
applied  in  the  image  domain.  Many  authors  propose  denoising 
techniques  by  assuming  a  Gaussian  noise  with  zero  mean  and 
variance  I  [10].  A  robust  estimate  of  the  noise  standard 
deviation  a^  is  obtained  in  the  finest  decomposition  scale  by 
measured  wavelet  coefficients  as: 


p{x\w,^)  = 


p{w\x,^)p{x\^) 


p{w\^) 

p{w\x,^)p{x\^) 


(9) 


Unknown  parameters  O  =  cr^,cr,^]  can  be  estimated 

using  EM  algorithm: 


log  p{x  I  w,  O)  = 

N  y 

Z S  I  [log  \x,,^)  +  log  p{x,  I O)  +  log  i]] 


(10) 


log;?(x|H’,6»)  = 
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(11) 


where  E[z.^  \  x\  is  defined  in  (19)-(20)  is  Gaussian  with 
0  mean  and  variances  al  and  al  -In  equation  (14)  unknown 
parameters  are  0=[x„  aP,  aP,  Oj}-  By  calculating  M-step 
the  estimates  (15)-(18)  are  obtained. 


<7= — [ — M4Z)({w,,  j,0<k<2')  (7) 

where  iv  denotes  a  wavelet  coefficient  and  J  denotes  the  finest 
level  of  wavelet  decomposition. 

B.  Expectation  maximization  algorithm  for  SAR  image 
denoising 

The  idea  is  to  estimate  unknown  parameters 
O  =  <jf  ,cr^]  using  EM  algorithm  and  wavelet 

decomposition  trees. 

In  this  case,  the  complete  dataset  consists  of  the  joint  events 
(xfo  Zk},  k=\,...,N  and  zt  takes  integer  values  in  the  interval  [1, 
J\,  and  it  denotes  the  mixture  from  which  x*  is  generated. 

P(.x)  =  YjPi.x\Zj)Pj 

^Pj=\\  p{x\  z.)dx  =  1  (8) 

>1  X 

Using  Bayesian  inference: 


Y,P[zi\x,9{t)] 
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'llalj(t)  +  llal.(t)\Y^P[z,\x,9(t)\ 
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wavelet  eoeffieient  and  predietion  error  is  eoded  using  image 
fuzzy  based  eoder. 


Y,PbAx,e{t)\[w,- 

Y^p[z,\x,e{^] 
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(15) 


IV.  EXPERIMENTAL  RESULTS 


P!{t  +  l)  =  ^f^P{z.\x„e(t)) 

^  k=l 


P{Zj\x„0{t))  = 


pix^\Zj,0{t))P.{t) 
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A.  Denoising  results 

The  X-SAR  image  256x256  with  2-looks,  shown  in  Fig. 2(a)  is 
transformed  with  dyadie  wavelet  transform  using  "symlets4" 
filter  bank.  The  image  is  denoised  using  proposed  method  as 
shown  in  Fig.  2(b).  3  levels  of  wavelet  deeomposition  were 
used. 


The  algorithm  repeats  until  ||0(?  -i- 1)  -  <I>(t)||  <  e  ■ 

The  dependenee  parent,  ehildren  in  the  dyadie  wavelet 
deeomposition  are  shown  in  Fig.  1.  Eaeh  parent  has  four 
ehildren  in  lower  deeomposition  level. 

The  eontext  seleetion  is  different  for  varianee  estimation 
al  and  crV  For  estimation  of  varianee  <j^  in  equation  (17), 

wavelet  coeffieients  that  surround  observed  eoeffieient  are 
eonsidered.  If  for  example  diagonal  wavelet  eoeffieient  on 
seeond  level  is  observed,  horizontal,  vertieal,  left  and  right 
neighbor  eoeffieients  are  eonsidered  and  the  same  teehnique  is 
used  for  eoeffieient  in  vertieal  and  horizontal  subband  at  the 
same  level.  The  observed  eoeffieient  is  predieted  using: 

V  =  +  6*1  (W^  + 

e^{Wu  +  Wa)  +  +  Wfj,,) 


Figure  1 .  Original  X-SAR  image. 


For  varianee  estimation  al  eorresponding  ehildren  wavelet 

eoeffieients  that  are  2  or  one  1  level  lower  as  observed 
eoeffieient  are  eonsidered.  In  Fig.  1  a  partieular  example  is 
shown,  where  eoeffieients  1  level  lower  of  observed  wavelet 
eoeffieient  are  eonsidered. 

III.  Compression  of  denoised  image 

Spaee  Frequeney  Quantization  (SFQ),  and  a  fuzzy  eontext 
based  modeling.  The  Spaee  Frequeney  Quantization  (SFQ)  is 
an  effeetive  ending  teehnique,  whieh  uses  a  zerotree  pruning  of 
wavelet  eoeffieients.  After  the  Space  Frequency  quantization  is 
applied  to  the  wavelet  coefficients  the  non-pruned  coefficients 
are  quantized  using  context  modeling.  The  context  based 
quantization  technique  uses  the  partitioning  proposed  by  the 
Trellis  Coded  Quantization  and  the  fuzzy  logic  to  predict  the 
next  TCQ  state  from  the  context.  Quantized  indices  are  entropy 
coded  using  an  arithmetic  coder.  The  probability  estimation  is 
based  on  the  observation  of  the  past  coded  bits. 

The  texture  parameters  obtained  at  the  despeckling  stage  are 
used  in  the  compression  process.  Texture  parameters  also 
defines  different  areas  of  the  image  (urban  scene,  agricultural 
scene,  etc.)  and  are  used  to  predict  an  observed  wavelet 
coefficient  from  the  past  wavelet  coefficients.  This  approach 
enabled  un-uniform  quantization  of  different  scenes  in  the 
single  image.  The  prediction  is  reduced  from  the  observed 


Figure  2.  Denoised  image  using  proposed  method. 


B.  Compression  results 

The  denoised  image  and  original  image  are  compressed 
using  compression  scheme  proposed  in  Section  III. 
Reconstructed  images  using  lossy  compression  method  are 
depicted  in  Fig.  3-6. 


Figure  3.  Original  image  compressed  to  1  bpp. 


Figure  4.  Original  image  compressed  to  0. 1  bpp. 


Figure  5.  Despeckled  image  compressed  to  1  bpp. 


Figure  6.  Despeckled  image  compressed  to  0. 1  bpp 


V.  CONCLUSION 

In  this  paper  despeckling  and  compression  algorithms 
based  on  wavelet  transform  are  presented.  Experiments  showed 
that  the  expectation  maximization  algorithm  can  be  used  in  the 
process  of  denoising  SAR  images.  Only  Gaussian  model  is 
considered.  It  is  well  known  that  wavelet  subbands  does  not 
have  Gaussian  distribution  thus  other  models  should  be  tested. 
Visual  impressions  of  compressed  despeckled  images  and 
original-speckled  images  at  higher  compression  ratios  are 
almost  the  same. 
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